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Abstract. We consider an incompressible fluid in a three-dimensional 
pipe, following the Navier-Stokes system with classical boundary conditions. 
We are interested in the following question: is there any optimal shape for 
the criterion "energy dissipated by the fluid"? Moreover, is the cylinder the 
optimal shape? We prove that there exists an optimal shape in a reasonable 
class of admissible domains, but the cylinder is not optimal. For that pur- 
pose, we explicit the first order optimality condition, thanks to adjoint state 
and we prove that it is impossible that the adjoint state be a solution of this 
over-determined system when the domain is the cylinder. At last, we show 
some numerical simulations for that problem. 
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1 Introduction 



The shape optimization problems in fluid mechanics are very important and 
gave rise to many works. Most often, these works have a numerical character 
due to the intrinsic difficulty of the Navier- Stokes equations. For a first 
bibliography on the topic, we refer e.g. to [7], [H, [H], [H] [16] . 

In this work, we are interested in one of the simplest problem: what 
shape must have a pipe in order to minimize the energy dissipated by a 
fluid? For us, a pipe (of "length" L) will be a three dimensional domain 0, 
contained in the strip {(^i, 3:2, X3) ,0 < X3 < L}. We will assume that the 
inlet E := dQ fl {x^ = 0} (where dQ denotes the boundary of and the 
outlet S := dil fl {xs = L} are two fixed identical discs and that the volume 
of Q is imposed. The unknown (or free) part of the boundary of Q will be 
denoted hy T [so dQ = E UT U S). 

In the pipe $1, we consider the flow of a viscous, incompressible fluid 
with a velocity u and a pressure p satisfying the Navier-Stokes system. We 
assume that the velocity profile uq at the inlet E is of parabolic type; on 
the lateral boundary F, we assume no-slip condition u = and we control 
the outlet by imposing an "outlet-pressure" condition on S. We will assume 
that the viscosity /j, is large enough in order that the solution of the system 
is unique (see p5|). The criterion that we want to minimize, with respect 
to the shape Q, is the energy dissipated by the fluid (or viscosity energy) 
defined by J{Q) := 2/i |e(u)pdx where e is the stretching tensor. 

We will first prove an existence Theorem. To obtain this result, we work 
in the class of admissible domains which satisfy an e-cone property (see 0], 
[9]). Then, we are interested in symmetry properties of the optimal domain. 
For the Stokes model, we are only able to prove that the optimum has one 
plane of symmetry. It is not completely clear to see whether the optimum 
should be axially symmetric. In a series of papers [2], [IS], G. Arumugam and 
O. Pironneau proved for a similar, but much simpler problem that one has to 
build riblets on the lateral boundary to reduce the drag. Nevertheless, it is a 
natural question to ask whether the cylinder should be the optimum for our 
problem. We will show that it is not the case. For that purpose, we explicit 
the first order optimality condition. This condition can be easily expressed 
in term of the adjoint state and gives an over-determined condition on the 
lateral boundary F. Then, we prove that it is impossible that the adjoint 
state be a solution of this over-determined system when the domain is the 
cylinder. 

This paper is organized as follows. At section [2l we state the shape 
optimization problem, we prove existence and symmetry. Section[3]is devoted 
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to the proof of the main Theorem. We give in section [4] some numerical 
results and concluding remarks. 

These results have been announced in the Note [lOj. 



2 The shape optimization problem 

Let us give the notations used in this paper. We consider a generic three 
dimensional domain contained in a compact set 



where Rq and L are two positive constants. We will denote by the 
boundary of fi. In the sequel, we will assume that the inlet E of 0, defined 
by E := dn n {x3 = 0} and the outlet S defined by S := dQ n {x^ = L} 
are two fixed identical discs of radius R < Rq centered on the xs axis. 
We will also assume that the volume of all the domains O is imposed, say 
|r2| = y = irR^L. We decompose the boundary of O as the disjoint union 
dQ = E U T U S and F, the lateral boundary is the main unknown or the 
shape we want to design. 

Let us now precise the state equation. We consider the fiow of a viscous 
incompressible fiuid into O. We denote by u = (ui,U2,U3) (letters in bold 
will correspond to vectors) its velocity and by p its pressure. As usual in 
fluid mechanics, we introduce e the stretching tensor defined by: 



We will consider the Navier-Stokes system (except for Theorem 12.41 where 
the Stokes system will be considered). As boundary conditions, we assume 
that the velocity profile uq at the inlet E = {x^ = 0} is of parabolic type; 
on the lateral boundary F, we assume adherence or no-slip condition u = 
and we control the outlet by imposing an "outlet-pressure" condition on 
5* = {x-s = L}. Therefore, the p.d.e. system satisfied by the velocity and 
the pressure is: 



D := {{xi,X2,X3) , xf + xl < RI ,0 < xs < L] 




(1) 



< 



-;uAu+ Vp + Vu • u = X € r2, 

div u = X G 

u = uo := (0,0,c(xf +xi - xG^; 

u = X G r 

—pn + 2fie{u) • n = h := {2ficxi, 2fj,cx2, —pi) x G 5. 
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where fi > denotes the viscosity of the fluid, n the exterior unit normal 
vector (on S we have n = (0,0, 1)). At last, the constant c which appears 
in the boundary condition on E and S is assumed to be negative. The sign 
of c can physically be explained. Indeed, in the case where O is a cylinder, 
the flow is driven by a Poiseuille law (simplified physical law derived from 
the Navier-Stokes system which describes a slow viscous incompressible flow 
through a constant circular section). Then , this constant c can be written 

c = — r^, where pi denotes the constant value of the pressure at the 

outlet S while po is the constant value of the pressure at the inlet E. 

This choice of the boundary condition ensures that the solution of ([I]) 
will be given by a parabolic profile when O is a cylinder. More precisely, if 
O is the cylinder of radius R and height L, the solution of ([T]) is explicitly 
given by: 



(2) 



u(xi, X2, X3) = (0, 0, c(xf + xl- R'^)) 

p{xi,X2,X3) = 4/ic(x3 - L)+pi. 



More generally, if is a regular domain, we have a classical existence and 
uniqueness result for such systems, see e.g. [3], [T9] . 



Theorem 2.1. Let us assume that uq belongs to the Soholev space {H^/'^{E))^ 
and h G {H^/'^{S)f. If the viscosity fj, is large enough, the problem (Op has 
a unique solution {u,p) G H^{Q,) x L'^{Q). 

The criterion we want to minimize is the energy dissipated by the fluid 
(or viscosity energy) defined by: 



(3) 



J{n) := 2/i / \e{u)\^dx, 
Jq 



where e is the stretching tensor : 



1 / dui duj 

2 \dxj dxi 



1<«J<3 



To make the statement precise, we also need to define the class of admissible 
domains or shapes. We will consider a first general class: 

(A) '= bounded and simply connected domain in : 

\n\ =V, iLo r\Jl = E, IiLr^Tl = S,] 

where Hq and XI^, denote respectively the planes {xs = 0} and {X3 = L}. 
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To prove an existence result, we need to restrict the class of admissible 
domains. It is a very classical feature in shape optimization, since these 
problems are often ill-posed, see [T|, [9]- We adopt here the choice made by 
D. Chenais in [4j which consists in assuming some kind of uniform regular- 
ity. More precisely, we will consider domains which satisfy an uniform cone 
condition, we say that these domains have the e-cone property, we refer to 
[1], [5] or [9] for the precise definition. So, we define the class 

(5) Oy := {0, G Ov ■ ^ has the e-cone property} 

Lemma 2.2. The class Oy is closed for the Hausdorff distance. 

Proof. We recall that the class of open sets with the e-cone property is closed 
for the Hausdorff convergence (see Theorem 2.4.10 in [9]). Moreover, the 
convergence also holds for characteristic functions, so the volume constraint 
is preserved. So, it remains just to prove that the properties defining the inlet 
E and the outlet S are preserved. Let {Qn)neN be a sequence of domains in 
Oy which converges, for the Hausdorff distance, to a domain Q. We want to 
prove that HqH^I = E and H/^ n = 5". The first inclusion Hq H C is 
just a consequence of the stability of inclusion for the Hausdorff convergence 
of compact sets. Let us prove the reverse inclusion: let xq G and n G N. 
Since fin has the e-cone property, there exists a unit vector such that the 
cone C(e,xo,i^ra) be contained in Qn- Up to a subsequence, one can assume 
that (^„) converges to some unit vector ^ and that the sequence of cones 
C(e,xo,Cn) converges (for the Hausdorff distance) to the cone C(e,xo,^). 
By stability with respect to inclusion, one has 



Therefore xq G 0,, and since xq G C Hq, the reverse inclusion is proved. 



where J is defined in |^ with u the velocity, solution of the Navier-Stokes 
problem (OP, and Oy is defined in has a solution. 




□ 



We are now in position to give our existence result. 
Theorem 2.3. The problem 
, , r min J(r2) 
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Proof. Let {^ln)n£N, be a minimizing sequence in Oy. Since the open sets 
fin are contained in a fixed compact set -D, there exists a subsequence, still 
denoted by Qn which converges (for the Hausdorff distance, but also for the 
other usual topologies) to some set Q. Moreover, according to Lemma [221 
O belongs to the class Oy. 

To prove the existence result, it remains to prove continuity (or lower- 
semi continuity) of the criterion J. For any n € N, we denote by Un and Pn 
the solution of the Navier-Stokes system ([T| on 0„. Due to the homogeneous 
Dirichlet boundary condition on the lateral boundary F, we can extend by 
zero Un and pn outside Qn- So we can consider that the functions are all 
defined on the box D and the integrals over Qn and over D will be the same. 
Let us first remark that (un) is uniformly bounded in H^{D). Indeed, the 
sequence \e{un)\'^dx = J'^|e(un)pdx is bounded by definition and the 
result follows using Korn's inequality on the set D together with a Poincare's 
inequality (see below proof of proposition I3.ip . 

Therefore, according to refiexivity of and the Rellich-Kondrachov's 
Theorem, there exists a vector u E [H^{D)]^ and a subsequence, still denoted 
Un such that : 

Un ^ u and Un — > u, Mq £ [1, 6[. 

It remains to prove that u is the velocity solution of the Navier-Stokes system 
on Let us write the variational formulation of ([1]). For any function w 
satisfying 

w e [H^{D)f : w = on S U F and divw = in D, 
and for all n G N, the function Un verifies : 



(7) 



/ (2^e(un) : e(w) + Vun • Un • w) dx = / h.Un • v^ds 
Jd Js 



Since we have weak convergence of Un, it comes 



/ e(un) : e(w)dx — 
Jd 



-oo 



e(u) : e(w)dx. 

D 



Let us now have a look to the trilinear term. We already know that Vun ^ 
Vu. Moreover, from Cauchy-Schwarz's inequality and Sobolev's embedding 
Theorem, we have: 



Un - u) • Vlf||r^2(£,)i3 < {Un,i-Ui)^dx / wfdx 




n 



l|2 II ||2 

"II [L4(D)]3 II W||[^4(£,)]3- 
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Then (un • w)„gN converges strongly in L'^{D) to u • w. Therefore, 



Vun • Un ■ wdz > / Vu • u • wdx. 

D n-»+oo 

Finally, weak convergence of Un in [H^{D)]'^ implies weak convergence of 
the trace in L'^{S) and the boundary term h.Un • wds in jT]) converges 
to /^h.u ■ wds. Therefore, u satisfies the variational formulation ^ (and 
also the boundary condition u = uq on £^ because every Un satisfies it). To 
conclude, it remains to prove that u is zero on the lateral boundary F. It is 
actually a consequence of the convergence in the sense of compacts of Vtn to 
f2, and the fact that O is Lipschitz and then stable in the sense of Keldys. 
We refer to Theorem 2.4.10 and Theorem 3.4.7 in [9]. □ 

We are now concerned with symmetry properties of the minimizer. When 
the state system is Stokes instead of Navier-Stokes the following result can 
be proved: 

Theorem 2.4. There exists a minimizer of the problem ^ (with the Stokes 
system as state equation) which has a plane of symmetry containing the ver- 
tical axis. 

Moreover, any minimizer of class has such a plane symmetry. 

Proof. Let O denotes (one of) the minimizer(s) of problem ^ and D the 
vertical axis xi = X2 = 0. Among every plane containing D, at least one, say 
Pq, cuts O in two sub-domains Qi and Q2 of same volume (volume equals to 
V/2). 

Let us now introduce the two quantities Ji and J2 defined by: 
Ji := 2fi / [e(u)pdx and J2 := 2fi / |e(upda;, 

so J{yi) = Ji + J2. Without loss of generality, one can assume Ji < J2. Let 
us now consider the new domain = ili \Ja{yti), where a denotes the plane 
symmetry with respect to Pq. We also introduce the functions (u,ij) defined 

by 

_ / if X G f^i _ r p(x) if X G f)i 

^ > 1 u(a(x)) ifxGa(Oi) ''^ > ~ \ p{<j{iC)) ifxGa(17i) 

It is clear that u G [H'^{Vt)f', p G L'^{VL) and divu = 0. Moreover 
2/x / |e(S)pda; = 4^ /" \e{u*)\^dx = 2Ji < J(0). 
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Now, it is well known that the solution of our Stokes problem can also be 
defined as the unique minimizer of the functional 

iPni^)) := 2fi [ |e(v)pdx 
Jn 

on the space 

V{n) := {v G H^{n) : divv = 0, V|^ = Uq and = 0}. 
Therefore, we have: 

< 2/u/f^|e(G)|2dx < J{n), 

this proves that Jl, which has the same volume as 0, and is symmetric with 
respect to Pq, is also a minimizer of J. 

Now, let us prove that if O is regular enough (actually but one can 
weaken as shown by the proof below), it must coincide with and there- 
fore is symmetric. Necessarily, we must have the equality in the chain of 
inequalities It proves, in particular, that u is the solution of the Stokes 
problem on Q. But since u coincides with u on Qi by definition, one can use 
the analyticity of the solution of the Stokes problem (see e.g. [T^) to claim 
that u = u on n 17. Now, if Q would not coincide with $7, we would have a 
part of the boundary of Q, say 7 included in Q. By assumption, Q being C^, 
the solution of the Stokes problem is continuous up to the boundary (see [8]) 
and therefore u should vanish on 7. By analyticity, it would imply that it 
vanishes identically: a contradiction with the boundary condition on E. □ 

As explained in the introduction, one can wonder whether the minimizer 
has more symmetry. In particular, could the cylinder be the minimizer? The 
following Theorem proves that it is not the case. It is the main result of this 
paper. The proof is absolutely not obvious and will be given at the next 
section. Let us remark that the following result also holds for the Stokes 
equation. The proof in the Stokes case follows the same lines and is a little 
bit simpler, see [17] for details. 

Theorem 2.5. The cylinder is not the solution of the shape optimization 
problem 

r minJ(r!) 

^ ' \ ne Ov, 

where J is defined in |^ with u the velocity, solution of the Navier-Stokes 
problem (OP, and Oy is defined in 
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3 Proof of the main theorem 

In all this section, Q will now denote the cylinder {x1+X2 < R^,0 < < L}. 
3.1 Computation of the shape derivative 

Let us consider a regular vector field V : M'^ ^ with compact support 
in the strip < x^ < L. For small t, we define fit = (/ + tY)Q, the image 
of Q by a perturbation of identity and f{t) := J{^t)- We recall that the 
shape derivative of J at O with respect to V is /'(O). We will denote it 
by dJ(r2; V). To compute it, we first need to compute the derivative of the 
state equation. We use here the classical results of shape derivative as in 
[9], [l3], [18]. The derivative of (u,p) is the solution of the following linear 
system: 

-/iAu' + Vu • u' + Vu -u + Vp =0 xen 



(10) 



div u' = 
' - 



u 



u 



= -T^(V-n) 
an 

p'n + 2fj,e{\i') • n = 



Now, we have (see [9], |18| ) 

(11) dJ{n,V)=4fi I e(u) : e(u')dx + 2/u / |e(u) 

Jn Jr 



X G 
xe E 

X G r 

X G S. 



(V • n)ds. 



It is more convenient to work with another expression of the shape deriva- 
tive. For that purpose, we need to introduce an adjoint state. 



Proposition 3.1. Let us consider (v, 
problem : 



(12) 



, solution of the following adjoint 
-2/iAu 



-^Av + Vu • V - Vv • u + Vq' = -2/iAu x G 17 

divv = X G 

v = xG-Eur 

—qn + 2/xe(v) • n + (u • n)v — 4;ue(u) • n = x G S. 



// the viscosity fi is large enough, then the problem H^) has a unique solution 
{'v,q)- Moreover, this solution belongs to C^{Q) x C^{Q). 

Proof. The existence and uniqueness of the solution is a standard application 
of Lax-Milgram's lemma. We introduce the Hilbert space 



V{n) := {u G H\n) : divu = 0}. 
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the bilinear form a and the Hnear form I defined by 



a(v, w) := / (2/ie(v) : e(w) + Vw • u • v + Vu • w • v) dx 

{£,w) := 4/i / e(u) : e(w)dx. 
Jn 

To prove ellipticity of the bihnear form a we use Korn's inequahty: 

l|Vv||[i2(n)p < Ci(||v|l[^2(f^)p + ||e(v)||[2.2(Q)]3). 
and a Poincare inequahty: 

(13) l|v|l[L2(f,)p < C2 / [e(v)|2dx. 

Jn 

These two inequahties yield (we also use the explicit expression of u given 
in ([2]) to estimate the integrals containing u): 

/ min(l,C2) , ,,„9 „A ,, ,,9 
a(v, v) > (^^-^i^ - |cl(i?2 + 2R)j ||v||f^,(^)]3. 

and a is elliptic as soon as /i > ^'^^^^^11^(1^(7^)^'^^'' • ^o"^) existence and unique- 
ness of the solution follow from a standard application of Lax-Milgram's 
lemma together with De Rham's lemma to recover the pressure. 

It remains to prove the regularity of the solution. The C°° regularity in 
O on the one-hand and on the smooth surfaces E, S and the interior of the 
lateral boundary T on the other hand is standard (cf. [8]). The only point 
which is not clear is the regularity on the circles ECiT and 5nr. To prove 
it, one can use the cylindrical symmetry which is proved later (without any 
regularity assumptions) in Theorem l3.3l This symmetry allows us to consider 
a two-dimensional problem in the rectangle (0, R) x (0, L) into the variables 
r = {xi + x^)^!"^ and X3. For that problem, one need to prove regularity at 
the corners (i?, 0) and (i?, V). For that purpose, one extends the solution by 
reflection around the line r = this leads to a partial differential equation 
in the rectangle (0, 2R) x (0, L) whose solution coincides with our solution 
in the first half of the rectangle. The regularity, up to the boundary, of 
the solution of this elliptic p.d.e. is standard and the result follows. □ 

Let us come back to the computation of the shape derivative. We prove 
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Proposition 3.2. With the previous notations, the shape derivative of the 
criterion J is given by 

(14) dJ{n,V) = 2^y" (^(u) : e(v) - |e(u)p) {Y.n)ds. 

Proof. Using Green's formula in pT]l . one gets 



dJ{n,V) = 4/i / e{u) : e{u)dx + 2fi [ \e{u)\^{Y .n)ds 
Jn Jr 

= -2/i / ((Au + Vdivu) • u')dx + 4/x / e(u) • n • u ds 
Jn Jan 

+2ii [ |e(u)[2(V-n)ds 
Jan 

Now, let us multiply the first equation of the adjoint problem (fT2ll by u' and 
integrate over one obtains 

-/i / Av • u'dx + Vq- u'dx + / (Vu)^ • v • u'dx 
Jn Jn Jn 

— / Vv • u • u'dx = —2fi / Au • u'dx. 
Jn Jn 

Using one integration by parts and the boundary conditions satisfied by u' 
and V, we get 

/ (2//e(u') • e(v) - Vv • u' • u + Vu' • u • v) dx 

Jn 

— / cr(v, q) ■ n ■ u'ds + / ((u • v)(u' • n) — (u • n)(u' • v)) ds 
Js Js 

— / cr(v, q) ■ n) ■ u'ds = —2^ / Au • u'dx. 
Jr Jn 

In the same way, if we multiply the first equation of the problem (fTO]l by v 
and integrate over we obtain 

—fi / Au' • vdx + / Vp' ■ vdx + / Vu' • u • vdx + / Vu • u' • vdx = 

Jn Jn Jn Jn 



and 



/ (2^e(u') • e(v) + Vu' • u • v - Vv • u' • u) dx 

Jn 

+ / {~(^{^'^P) • n • V + (u • v)(u' • n)) ds = 0. 
Js 
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Coming back to the shape derivative expression 

dJ{n, V) = -2/i / ((Au + Vdivu) • u')dx + 4// / e(u) • n • u'ds 
Jn Jan 

+2fi [ |e(u)|2(V-n)ds 
Jan 

= A + 4fi e(u) • n • u'ds + 2/i / |e(u)p(V • n)ds, 
Jan Jan 

where we set A := — 2/i / ((Au + Vdivu) • u')dx. Using the previous iden- 

Jn 

tities, we get for A 

A = (gn - 2;ue(v) • n) • u'ds - / (u • n)(v • u')ds. 

Jrus Js 

Therefore, according to (fT2]l 



dJ{n, V) = /" {qn - 2/ie(v) • n) • u'ds - / (u • n)(v u')ds 

Jrus Js 

+4:fj. / e(u) • n • u'ds + 2/i / |e(u)p(V.n)ds 
Jsur Jr 

- j (gn - 2/ie(v) • n + 4/xe(u) • n) • u'ds + 2/i j |e(u)p(V.n)ds 
j (^{qn - 2/ie(v) • n + 4/ie(u) • n) • |^ - 2/i|e(u)|2^ (V • n)ds 



To get the (more symmetric) expression given in (|14|) , one can use the follow- 
ing elementary properties. Since u (and v) is divergence-free and vanishes 
on r, we have on this boundary: 

. e(u).n.f^ = |e(u)|2. 
. (e(v).n).|i=e(u):e(v). 
Proposition 13.21 follows. □ 

3.2 Analysis of the PDE M l 

We will prove the following symmetry result for the solution of the adjoint 
system. It shows that the solution has the same symmetry as the cylinder. 
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Lemma 3.3. 

With the same assumptions on ^ as in Proposition \3.1\ there exist {w,W3) G 
[iJ^((0,fi)x(0,L))]2 andqe L^{{0, R)x{0, L)) such that, for any {xi,X2,X3) G 

n 

(i) Vi{xi,X2,X3) = Xiw{r,X3), for i G {1,2} ; 

(ii) V3{xi,X2,X3) = W3{r,X3) ; 

(iii) q{xi,X2,X3) = q{r,X3). 
where r = (xf + Xg)""^/^. 

Proof. Let us introduce the differential operator 2g defined by 

d d 

= Xl^ X2T, — • 

0x2 0x1 

£0 corresponds actually to the differentiation with respect to the polar angle 
0. Let us set 

(15) Vi = ^e{v^), ViG {1,2,3} andg = £e(g)- 

By applying the operator £e to the equation (fT2]l we get the following system 
(where we have used the explicit expression of the solution u given in ([2])) 

(16) ^ 

-fi Avi + 2cxiV3 - 2CX2V3 - c(xl + xo - -R^)7r^ + — = x G S7 

0x3 oxj 0x2 
dv dQ dQ 

-fi Av2 + 2CX2V3 + 2cxiV3 - c(x? + Xn - R^)t[ 1- -7[ 1- -7[ — = X G O 

^ ^ 0x3 0x2 OXl 

-nAiJ3-c(xl + xl-R^)^ + -^ = xGfi 

0x3 dx3 

^ + ^ + ^-^ + ^ = X G 

dxi dx2 8x3 8x2 8x1 

Vi=V2 = V3 = Q xG-Eur 

/i f + — ) - //tt-^ + c(x? + Xo - R^)vi = -A11CX2 X G S, 

\8x3 0x1 I 0x2 

( dv2 8V3\ dv3 2 , 2 ry2^^ . 

fi h ^ — + /i^ h c(xi + X2 - -K )f2 = 4;ucxi X G 6, 

Voxa ax2/ axi 

2/u|^ + c(x? + x| - i?2)w3 = g X G 5", 

0x3 

Let us now introduce the following new functions 
• zi = vi + V2 ; 
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• Z2 = V2-Vi ; 

• Z3 = V3. 

According to system (fT2]l . the system (fTHjl rewrites in term of zi, Z2, Z3 
-fj,Azi + 2cxiZ3-c{xl + xl- R^)- h TT-^ = xen 



(17) 



0Z2 



dxj^ 
dq 



-fiAz2 + 2CX2Z3 - c{x\ + + — ^ = xeO 



3X3 9X2 



A / 2 , 2 ry2^dz3 dq 

9X3 

Zl= Z2 = Z3 = 
a£2 _^ dz3_ 

^8x3 dx2 
2/i|^ + c{xl + xl- R^)z3 

0x3 



+ Z2c{xl + X2 







X G r2 
X G 

X G -Bur 

xG 5, 
xG 5, 
xG 5, 



This adjoint problem has a unique solution if ^ is large enough (see propo- 
sition [3]T]) , therefore 

zi = Z2 = V3 = q = 0■ 
The fact that ^3 = 2e{v3) and q = £,e{Q) vanish proves points ii and iii of the 
Lemma. Now let us precise the properties of functions vi,V2- It has been 
proved that 2^0 (vi) = —V2 and £^0(^2) = vi. Therefore, applying once more 
the operator 2g yields £9 ° 2e{vi) + vi = 0. This implies that there exist 
two functions a and (3 in the space H^{{0,R) x (0, L)), such that 

vi = xia{r, X3) + X2/3(r, X3). 

Moreover, since ileivi) = —V2, we get 

V2 = -xi/3(r, X3) + X2a(r, X3). 

To finish the proof, it remains to check that the function (3 is identically 
zero. For that purpose, let us write down the partial differential equation 
satisfied by j3. Prom the two first equations of system (fT2]l and the boundary 
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condition, we can prove that /? satisfies the following system 
(18) 

-4 a^l + + ^ j - - = -3) ^ (0, R) X (0, L) 

< /?(r,0)=/?(i?,X3) = |^(0,X3)=0 (r,X3) G (0,i?) X (0,L) 

+ c(r2 - i?2)^ = (r,X3) G (0,i?) X {L} 

It remains to prove that the zero function is the unique solution of the 
previous system. Multiplying the equation by /3 and integrating on the 
rectangle in polar coordinates gives, using the boundary conditions 

+/i / /32(0,:E3)dX3 + - / ir'' -R^)f{T,L)TAT. 

Jo ^ Jo 

Since c < and r < /Z, we get |f = in {0,R) x (0,L) and (3'^{0,X3) = 
for any X3 € (0, L). Then /3 = which gives the desired result. □ 



3.3 The optimality condition 

We argue by contradiction. Let us assume that the cylinder 17 is optimal 
for the criterion J. We first write down the first order optimality condition. 
From the explicit expression ([2)) of u, we have 

/ cxi \ 
e(u) = CX2 . 

\ CXl CX2 / 

Therefore 

\e{u)f = 2c\xl + xl), 

and |e(u)p = 2c^B? is constant on F. 

Now the first order optimality condition ensures the existence of a La- 
grange multiplier A G M, such that dJ(r2, V) = AdVol (17, V) for any vector 
field V. Due to the expression of the shape derivatives of J and the volume, 
it writes 

2^1 J (e(u) : e(v) - |e(u)p) (V.n)ds = A y" (V • n)ds. 
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This implies that e(u) : e(v) is constant on T. Now, from the expression of 
e(u) on r, we deduce 

c f dvs dvs dvi dv2\ 
c / \ cR dvs 

= 7^ TT h X2- 



' dxi dx2 J 2 5n |p' 
because vi\-^ = ^2[p = 0. Therefore the optimahty condition writes 

(19) 3^ G M : ^ = e on r. 

on 

Now, we give another useful Lemma 

Lemma 3.4. // the cylinder is optimal and using the notations of Lemma 
Ig.gj. we have 

dq_ ^dq ^ ^ 

Proof. Let us write the adjoint problem (fT2]l in term of the functions w, ws 

et q. We get 

(20) 



+ -ir + irr + -ir + ^c^^s - c{r' - R)^ = o in o 

\ ar^ r or oxi ) r or 0x3 
^/ (9^W3 ^ ldw3 ^ d^W3 \ ^ 1 dq 2 ^2n^'^'3 _ ^ 

dw dW3 ^ 
2w + r— + — i = m 

Or ax 3 

u;(r, 0) = W3{r, 0) = u;(i?, X3) = W3{R, X3) = 
^ dw ldw3\ 2 D2^ 



/X — + — + c(r^ - R'')w = 4//C on 5 
\t)X3 r or J 

o , / 2 E>2x ~ c 

ZyU— h c(r — R )W3 = q on b. 

0x3 

Since w\, = tt'si, = 0, we have S^, = = and 

= 0. In particular, from the divergence-free condition, we obtain 

^3 \{r=R} 

dw Q 

Now, let us differentiate the divergence-free condition with respect to r, 
we get 

V(.,X3)E(0,i?)x(0,L), 3|^+rf| + |^ = 0. 

Or Or'' Orox3 
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Now, ^ = ^ (it is the optimality condition (fT9l) ) ; therefore, we have 

a I = 0. Combinine; this last result with = 0, it comes 

9^ If; 



0. 



We let r going to R in the first equation of problem ([20]) and we use the 
previous identities to get 

□ 

3.4 An auxiliary function 

Using notation of Lemma 13.31 we introduce now two new functions 

• Wo : [0, R] X [0, L] — > R 

(r, X3) I — > / w{r,z)dz 
Jo 

• tp: [0, i2] X [0, L] — > R 

(g(r, X3) — 2cr^wo(r, X3)) d^rdr 

We will also denote by the horizontal section of the cylinder {x € : 3:3 = z}. 
The following lemma is the key point of the proof. 

Lemma 3.5. The function ip is affine. 

Proof. The couple (v, q) satisfies the following p.d.e. 

-^Av + Vg + Vu • V - Vv • u = -2/iAu. 

Let us compute the divergence of both sides of the previous equality. Using 
the expression of u in the cylinder we obtain that (v, q) verifies 

21 ^q + \cv^ + 2c xi-^ + - 2c xi-^ + = 0. 

V dx\ 6x2 ) \ oxz 0x3 J 

Let us integrate this equation on a slide 

CO := {(xi, X2, X3) G ri; z- < X3 < z+} 
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(we will denote by e the inlet of uj and s its outlet). We get 

f . , , f f dv3 dvs \ f dvi dv2 \ , 

/ Aq + 4cf3dx + 2c / xi^r h X2t; — - 2c xi- h X2t; — dx = 0. 

Jcj Juj V dxi dx2j V ^^3 9X3 7 

Now, from Green's formula, we have 

/ ^I'^^^dx = / xi^snids— / fsdx = / xi^snids— ( fadx = — ( i^adx 

Ju) (yXi Jg^ 7a; J dujClF Ju Ju 



Therefore 



in the same way / X2Tr^dx = — / v^dx . 

OX2 



4c / v^dx + 2c I [ ^1^^ + ^2^^ ) dx = 0, 



so 



(22) / Aqdx = 2c I ( + X2^ 1 dx. 



3X3 C^2;3 ^ 

Let us consider the left-hand side of ((221) . Prom Lemma [3]4] it comes 

(23) / Agdx= / p-ds. 

Now, let us consider the right-hand side of ([22]) . Integrating by parts yields 

• / xi-^-^dx = / xivin^ds = f xifin3ds. 

• / a;2-^^dx = / X2V2n3ds = [ X2V2nsds. 

Jul C'X3 Jq^ JeUs 

Combining this result with ([23]l gives 

(24) ~ ^'^(^I'^i + ^2^'2)^ j ~ 2c(xit>i + X2U2) j ds, 

what can also be rewritten for any (z_,z_|_) G (0, L)^ : 
(25) 

-(r, z_) — 2cr'^w{r, z_) I rdr = / ( — — (r, z+) — 2cr'^w{r, z^) I rdr. 



dx3 J Jq \dx3 
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Now, since ip{z) = 27r / (g(r, z) — 2cr'^wo{r, z)) rdr, we have by differenti- 
Jo 

ating, for all z in [0,L], 

i^'iz) =271 I ( — 2cr^— — ) rdr = Itt ( 2cr'^w ] rdr. 

Jo \ox3 8x3 J Jo \dxs J 

Now, identity (f25]l proves that ■i/'' is a constant function which gives the 
desired result. □ 

We are now in position to precise the value of the constant ^ appearing 
in the first order optimality condition (fT9]) . For that purpose, we use the 
symmetry result given in Lemma 13.31 together with equation (l20l) . In this 
equation, let us integrate between X3 = and X3 = z G (0,L). Since 
wsir, 0) = 0, we get for any (r, z) £ [0, R] x [0, L] : 

dvjQ 

2wo[r, z) + r— — (r, z) + Ws[r, z) = 0. 
or 

Let us differentiate this last relation with respect to r. This yields 

. , dwo d'^wn dws 

26 3—^ + — ^ + ^ = 0. 

or or'' or 

Now, in ((20]) . we differentiate the divergence equation with respect to r, and 
we make r ^ R. We obtain 



dr |r dr"^ \ 



r 



Letting r going to R in (I26p and interverting limit and integral gives, using 
the previous equality 

^ =0. 

on |r 

So we conclude that = and the optimality condition rewrites 
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3.5 End of the proof 

Let us use the function ■0 defined above. We can rewrite it as 



0(z) 



/ [q — 2cr^ti;o) d^rdr = 2tt I {q{r, z) — 2cr'^wo{r, z)) rdr, 
Jt^ Jo 



where denotes the horizontal section of the cyhnder of cote z. We proved 
in Lemma 13.51 that is affine, therefore its derivative ^p' is constant, say 
^l^'{z) = a. The contradiction will come from the computation of this con- 
stant on the inlet E and the outlet S. We will see that we obtain two different 
values. Let us denote by A2 the two-dimensional Laplacian (with respect to 
the variables xi and X2)- 

Computation of the constant on the outlet S of the cylinder. First of all, let 
us remark that if we differentiate with respect to xi the boundary condition 
on S satisfied by the function vi, we get 

(28) fi + /U^-^ + 2cxivi + c{xl + xl- R^)^^ = 4^c, on S. 

OX1OX3 oxf 0x1 

In the same way, if we differentiate with respect to X2 the boundary condition 
on S satisfied by the function V2, we get 



Summing the two relations (1281) and (l29]) and using the divergence-free con- 
dition yields 

-J-^-Y + + 2c{x\V\ + X2V2) - c{x\ + X2 - i?^)^^ = 8/UC on S. 

Now, according to (fT2l) . ^3 satisfies 



(30) /XA2.3 = 8/xc - - c{x, + ^2 - ^ )^ + 
Combining together the two previous equations, it comes 

(31) - 2ij^ - 2c{xl + xl-R^)^ + p- + 2c{xivi + X2V2) = on 5. 

ox^ 0x3 axs 

Now, we integrate on S the equation (I30p . we have 



d'^vs 5^3 2 , 2 9q . I 

-fiA2V3 - l^-g^ - ^(^1 +^2- R )c + ^ ) ds = -8/ic / ds. 
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In the Proposition 13. 11 we have seen that V3 is up to the boundary. Taking 
into account the boundary condition on 5, we have 

A2Vsds = [ ^^dcr = . 
Jsnr on 

So, the integration gives 

-c [{xl + xl- R')^d.s + / p-ds = -Spic^R\ 
dxl Js 0x3 Jsdx3 

Using ([3T]) . we can deduce that 

1 f dq f 

- / ^ — ds - c {xivi + X2V2)ds = -SficirR^. 
^ Js Js 

According to Lemma [331 one can write 

xivi + X2V2 = {xf + xl)w [xl + xl) ^^"^ , X3 
Therefore 

(32) a = tp'{L) = -IQficTTR^ 

Computation of the constant on the inlet E of the cylinder. Let us first 
(9f 3 

remark that — — = (just use the divergence-free condition extended to 

0x3 |g 

E and the fact that vi\^ = V2\^ = 0). Let us now integrate the p.d.e. (fT2]l 
satisfied by ^3. We have, using |||| = 0, 

— /i / Awsds + [ -^^ds = — 8/ic / ds. 
Je Je 0x3 Je 

Taking into account the condition (f27ll we get 



-/X [ Av3ds = -fi [ A2V3ds - /i / ^r|ds 
Je Je Je dxi 



-'iE.r^^''^nE\d^.^d^.^^' 

( dvi dv2 \ 
/ ^^i — '^i + ^ — do- = 0. 
lEnT \OX3 0x3 J 



Then, it follows 



(33) / ^ds = -SucttR^. 

E 0x3 
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At last, since = V2\^ = 0, we have 

V''(0) = 2tt [ f P(r,0) - 2cr'^w{r,0)] rdr = [ -p-ds. 
According to ([33]l we have 

(34) a = i)'{<d) = -SficnR^. 

which is clearly a contradiction with ((32]) since c < 0. This finishes the proof 
of Theorem I2.5[ 

4 Some numerical results 

In this section are presented some numerical computations. It gives a con- 
firmation that the cylinder is not an optimal shape for the problem of min- 
imizing the dissipated energy. In particular, we are able to exhibit better 
shapes for this criterion. All these computations have been realized with the 
software Comsol. 

For any bounded, simply connected domain Q in or and any real 
numbers fi, b {b will be fixed in all the algorithm), let us define the augmented 
Lagrangian of our problem ^ by 

c{n,n) = + fi{\n\ -v) + ^ - y)^. 

Since Theorem 12.51 ensures that the cylinder is not optimal for the cri- 
terion J, the question of finding a better shape in the class of admissible 
domains Oy is natural. The numerical difficulties in such a work, are the 
non linear character of the state equation and the need to take into account 
the volume constraint. 

For that reason, we decompose the work in two steps. First, is considered 
a gradient type algorithm in two dimensions which allows us to reduce the 
criterion J. Then, we work in a three dimensional class of domains with 
constant volume V and cylindrical symmetry. In this class, we are able to 
find a shape (probably not optimal) which is better than the cylinder, see 
section 14.21 

4.1 A numerical algorithm in 2D 

We denote by the cylinder with inlet E, outlet 5, and measure V. ^Iq is 
our initial guess for the gradient type algorithm we consider. We deform 
by using the following method: 
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1. 



We fix /Uo G M, r > and e > 0. 



2. Iteration m. At the previous iteration, fj,m and 0^ have been com- 
puted. We define Om+i := (I + £mdm)(^m), where I denotes the 
identity operator, Em is a real number (step of the gradient method) 
which is determined through a classical ID optimization method and 
dm is a vector field of M^, solution of the p.d.e. 



where denotes the lateral boundary of O^, i.e. := d^mXiEUS). 
The solution of this p.d.e. gives a descent direction for the criterion J 
(see for instance [l], [6]). 

Then, the Lagrange multiplier fim is actualized by setting 



3. We stop the algorithm when (/im)m>o has converged and the derivative 
of the Lagrangian is small enough. 

The Figure [1] shows the geometry we obtain. The criterion has decreased 
about 1.1 % from the initial configuration (a rectangle here). 




X G 

X G ^US 

X G Tj^; 




Figure 1: Final 2-D shape obtained by the gradient algorithm 



23 



4.2 Some 3D computations 



In this section, we create a family of 2D shapes, constructed with cubic sphne 
curves which look like the presumed optimum obtained in figure [H Then, 
we obtain a family of 3D domains of volume V, by revolving the previous 
2D shapes around the (Oxs) axis. We introduce a small parameter e in the 
control points of the cubic splines and we evaluate for each value of e the 
criterion J. The value e = corresponds to the cylinder. Let us respectively 
denote by J(e) and J{0,q) the values of the criterion J evaluated at the 
domain corresponding to value e of the parameter and at the cylinder. Figure 
[2] is the plot of function e i-^ 100. above, and Figure [3] represents a 

better shape than the cylinder for the criterion J which is obtained with a 
value of the parameter e ~ 0.001. It shows that this simple method provides 
a 3D (axially symmetric) shape which is slightly better than the cylinder. 




0.001 0,002 0,003 0,004 0,005 0,006 0,007 0.008 0,003 0,01 



Figure 2: The cost function (which slightly decreases before increasing) 
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